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ABSTRACT 


A  discrete  model  suitable  for  the  analysis  of  polycrystalline  aggregate 
response  under  macroscopically  uniform,  quasi-static  loading  is  developed, 
with  particular  emphasis  on  the  characteristics  of  subsequent  yield  surfaces 
in  stress  space.  Internal  stress  and  deformation  states  are  determined  from 
approximating,  piecewise  linear  infinitesimal  displacement  fields  within 
crystal  grains,  based  upon  broadly  defined  constitutive  behavior  which  per¬ 
mits  inclusion  of  cubic  or  hexagonal  crystal  anisotropy  and  relatively  general 
hardening  laws  over  crystallographic  slip  systems.  Appropriate  aggregate 
matrices  are  established  as  symmetric,  positive-definite,  and  internal  fields 
corresponding  to  the  solution  of  the  discrete  model  are  proved  to  be  unique. 
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1.  INTRODUCTION 
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The  first  satisfactory  theory  for  predicting  the  plastic  deformation  of 
polycrystalline  aggregates  from  phenomenological  laws  of  single  crystal  behavior 
was  advanced  by  Sir.  G.  I.  Taylor  [1,  2].  In  Taylor's  now  classic  work,  the 
simplest  possible  kinematic  model  was  adopted  consistent  with  the  concept  of 
a  deformed  continuum  —  uniform  strain  throughout  the  crystal  grains.  This 
theory  was  generalized  by  Bishop  and  Hill  [3-5]  to  enable  the  approximate 
calculation  of  macroscopic  yield  surfaces  of  pronounced  yielding  (neglecting 
elastic  behavior)  and  modified  by  Lin  [6]  to  incorporate  elastic  strains. 
Subsequent  studies  were  made  by  Payne,  et.  al.  [7,  8].  Other  theories  and 
models  of  interacting  crystals,  all  utilizing  isotropic  elastic  field  solutions 
in  one  form  or  another,  have  been  proposed  and/or  Investigated  by  Krdner  [9], 
Eudlansky  and  Wu  [10],  Hutchinson  [11,  12],  Hill  [13],  and  Lin  and  his 
associates  [14-23],  with  the  models  of  Lin,  et.  al.,  most  nearly  satisfying  all 
equilibrium  and  kinematic  conditions  in  numerical  evaluations. 

In  the  present  paper,  a  new  discrete  aggregate  model  suitable  for  predicting 
the  response  of  thin-walled,  polycrystalline  tubes  is  presented  which  incorpo¬ 
rates  certain  features  of  previous  models  but  is  more  general  in  several 
Important  respects.  Moreover,  the  model  is  closely  related  to  theoretical 
characteristics  of  crystal  and  aggregate  behavior  established  by  Hill  [24,  25]. 
The  paper  is  organized  as  follows.  In  Section  2,  a  continuum  mathematical 
model  of  an  aggregate  of  identical,  polycrystalline  "unit  cubes"  is  introduced 
and  macroscopically  uniform  fields  are  defined.  The  aggregate  virtual  work 
equation  relating  microscopic  and  macroscopic  tensor  variables  then  follows 
as  a  direct  consequence.  Several  different  internal  stress  and  infinitesimal 


1 

i 


scraln  fields  are  presented  in  Section  3,  and  an  inequality  is  established  which 
leads  to  a  proof  in  Section  4  of  a  macroscopic  Bauschinger  effect  for  the  poly¬ 
crystalline  aggregate.  In  Section  5  a  further  specification  of  the  continuum 
model  is  given  appropriate  to  the  formulation  of  aggregate  boundary  value 
problems.  Sections  6  and  7  contain  the  general  analysis  of  the  discretized 
model  for  prescribed  aggregate  macrostrain,  based  upon  approximating,  piecewise 
linear  Infinitesimal  displacement  fields  within  crystal  grains.  Kinematic 
conditions  are  met  identically;  equilibrium  between  grains  is  satisfied  in  an 
average  sense  at  each  crystallite  node;  and  both  anisotropic  crystal  elasticity 
and  fairly  general  crystal  hardening  matrices  are  Included  as  constitutive 
behavior.  Aggregate  matrices  to  be  Inverted  (or  decomposed)  are  established  as 
symmetric,  positive-definite;  and  in  Section  8  a  strict  uniqueness  proof  is  pre¬ 
sented  for  both  Incremental  plastic  shears  in  crystallographic  slip  systems  and 
internal  stress-  and  strain-increment  fields  (as  determined  from  the  discrete 
modal).  Lastly,  in  Section  9  the  necessary  steps  for  calculation  of  subsequent 
yield  surfaces  era  given  and  a  suggested  model  for  quantitative  studies  is 
briefly  discussed. 


2.  PRINCIPLE  OF  VIRTUAL  WORK  WITHIN  THE  POLYCRYSTALLINE  AGGREGATE 


Consider  an  arbitrary  volume  V  of  surface  S  within  e  polycrystalline  metal 


specimen.  We  denote  any  statically  admissible  stress  field  in  V,  corresponding 

•  * 

to  a  system  of  salf-squllib rating  tractions  £  on  S,  by  £  (with  tensor  com¬ 
ponents  c*  )  and  any  continuous,  plscevlsa  differentiable  lnflnltaslsal  dls- 


plac 
yields 


•ij 

t  field  by  6u°.  A  straight forward  application  of  Gauss'  theorem  then 
% 


wherein 


/(;*•  6£°)'!V  -  / (T**  6u°)dS 

y  ^  ^  g  ^  *V 


„  iMrf'ieurrirti-  T-n-f  •  *-» 

V,r«^:^VW «' 

3 

(2.1) 


tro  «T  ,  o 
6£  •  P  6u 


(2.2) 


locally  vlchln  each  crystal  grain  of  V.  (See  the  Appendix  for  definitions.) 
To  obtain  an  equivalent  expression  for  the  right-hand  side  of  (2.1)  in  terms  of 
macroscopic  stress  and  strain  over  the  smallest  possible  volume  sample,  we 
introduce  the  mathematical  model  of  an  aggregate  of  identical  "unit  cubes"  and 
define  macroscoplcally  uniform  fields 


6u°(A+) 


6u°(A.)  +  c° 
\  1 


X*  + 

i  <*;>  - 


i*  - 

J  <V 


(2.3) 


(2.4) 


(respectively  kinematically  and  statically  admissible).  A^  denotes  the  unit 

cube  face  corresponding  to  the  positive  coordinate  axis  x.  and  c°  is  a  constant 

l  \ 

vector  independent  of  position  over  A+.  Macroscopic  stress  and  incremental 
strain  are  evaluated  in  a  natural  way  as 


*  .  f 

1J  At 

Ti  “i 

(no  sum 

station) 

(2.5) 

i 

4c°  •  — 
ij  2 

*.  „  J. 

(ci  +  ci 

)  . 

(2.6) 

s 


A 


Upon  substitution  of  (2.3),  (2. A)  and  (2.5)  into  (2.1),  we  find 


(2.7) 


From  considerations  of  moment  equilibrium 


A1 


/  I,  iA 


j 


(2.8) 


Hence,  the  macroscopic  stress  tensor  is  symmetric  and,  from  (2.6)  and  (2.7), 
(2.1)  can  be  written 


/  (?*  •  6£°)dV  -  o*.  6e°  (2.9) 

ith  V  now  representing  the  unit  volume  of  the  polycrystalline  cube.  Further¬ 
more,  from  (2.2),  (2.3)  and  (2.6), 

6e°  -  /  (65°)dV  (2.10) 

%  y  'V 

and,  from  (2. A)  and  (2.5)  and  force  equilibrium  over  any  interior  plane  area 
normal  to  a  coordinate  axis, 

a*  «  /  ?*  dV  .  (2.11) 

%  y  % 

Equations  (2.10)  and  (2.11)  are  of  course  equivalent  to  (2.6)  and  (2.5),  and 
aquation  (2.9)  (togather  with  these  definitions)  is  the  well  known  virtual 
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work  equation  for  the  polycrystalline  aggregate.  We  merely  remark  here  that 
by  mathematically  defining  states  of  homogeneous  macrostrain  and  macrostress 
(2.3,  2.4)  the  relationship  (2.9)  follows  immediately  and  the  various  additional 
arguments  of  Bishop  and  Hill  [3],  Kocks  [26],  and  Hill  [25]  ire  unnecessary. 
Bishop  and  Hill's  criterion,  which  can  be  written 


for  arbitrary  i,  j,  k,  is  in  fact  distinctly  different  from  (2.3,  2.4),  and 
neither  condition  implies  the  other.  In  addition,  although  the  macroscopic 
stress  tensor  is  symmetric,  equations  (2.3)  and  (2.4)  do  not  preclude  the 
existence  of  small  macroscopic  couple  stresses  (depending  upon  the  distribution 
of  crystal  orientations  within  the  cube).  These  are  determined  as 


where  r  is  the  position  vector  to  a  point  on  the  face  A 


(2.13) 


3.  SOME  GENERAL  INEQUALITIES  AND  INTERNAL  FIELDS 

Let  £  denote  the  local  stress  field  within  a  crystal  grain  determined  by 

assuming  that  the  aggregate  response  to  macrostress  o  is  wholly  elastic.  As  in 

(e ) 

[25],  we  take  ;  to  be  expressible  in  terms  of  o  and  a  tensor  (matrix)  function 

%  'V 

V  of  position  within  V  (i.e. ,  the  influence  of  elastic  inhomogeneity  on  the 
stress  field) : 


<?  o 


I 


5 


(e) 


<\i 


W  o 


6 

(3.1) 


Then,  from  (2.11) 


a  -  /  ?(e)dV  (3.2) 

^  y 


so  chat  /  Y  dV  -  I  (the  identity  matrix).  We  further  define  a  microstress  field 
S  (a}  v 

r  v  *  due  to  internal  slip  and  self-stressing, 

% 


;S(o) 


(3.3) 


with  (2.11  and  3.2) 


/  ?S(o)dV  -  0  ,  (3.4) 

\  % 

and  introduce  kinematically  admissible,  infinitesimal  displacement  and  corres¬ 
ponding  strain  fields  6uS^°\  such  that 

'V  'V 

55  -  C  5?(e)  +  55S(o)  .  (3.5) 

^  'V  r\j  i\, 

is  the  positive-definite,  elastic  compliance  matrix  of  an  individual  grain 
and  5,  5  are  the  actual  local  stress  and  strain  fields.  We  also  have 

Oi  r\j 


65  -  C  6c  +  65p 

r\j  'V  %  1 


(3.6) 
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g 

where  C  6c,  ■  6£  Is  Che  micro-elastic  strain  increment.  The  local  micro-plastic 

f\,  %  'Xi 

strain  Increment  6£P  is  determined  from  the  Incremental  plastic  shears  &Yk  on 

'V, 

the  N  slip  systems  of  the  crystal  through  the  transformation 


6£P  -  NT6y  ,  0.7) 

*\j  “V 

with  the  resolved  shear  stress  in  these  slip  systems  evaluated  as 


t  -  N  ;  .  (3.8) 

r\j  t\j  i\j 

The  transformation  matrix  N  is  defined  in  terms  of  its  kth  tow  vector  in  the 
Appendix.  (Opposite  eeneee  of  ellp  in  the  ease  crystallographic  slip  system 
are  denoted  by  different  k'e  eo  that  <Yk  is  always  non-negative). 

We  now  Introduce  the  following  scaler  averages  over  the  unit  cube: 


d  ■  /  («t  •  UP)dV 
P  y  1  v 


*  -  /  («C(,)*  6CS(0))dV 

P  y  %  \ 


/  (6;  •  C  6;)dV 

y  \  ^ 


(3.9) 


/  C  6;w)dV  . 

y  A/  'V  "V 


(e), 


(Note  that  (2.9)  applies  separately  only  to  w_  end  w_  since  6CP  end  C  5c  ere 

p  e  v  \  \ 

not  derivable  from  displacement  fields,  hence  are  not  separately  kinematically 


a 


t 


adalaalbla.)  Proa  (2.9),  (3. 4-3. 6),  and  (3.9) 


d  ♦  u  •  6o  •  6c  •  w  ♦  w 
P  •  x  -v  P  • 


u  -  w  •  /  (6CS<0>*  C  «tS(o))dV  >  0. 

•  •  y  \  V 


Hanca,  combining  (3.10,  3.11)  irlch  (2.9)  and  (3. 3-3.6), 


d„  -  *  •  /  (6cS<0> •  UP)dV  <  0 

9  P  y  %  % 


or  (ualng  3.7  and  3.8) 


8 

(3.10) 

(3.11) 

(3.12) 


/  (6tS(o>  •  6v  )dV  <  0  (3.13) 

y  \  * 

from  which  va  coo  dud  a  chat  cha  lacraaancal  ahaar  acraaaaa  It**3*  duo  co  alip 
and  aalf-acraaalng  oppoaa  cha  lacraaancal  plaadc  ahaar  a  4tk  la  a  aajorlcp 
(If  aoc  all)  of  cha  acclva  allp  apacaaa  of  cha  aggrogaca.  Thla  raaulc  will  ba 
cal  lad  upon  la  dacaralalag  canals  gaoaral  characcarladca  of  aubaaquanc  plaid 
aurfacaa  la  Sacdoa  4. 

la  alallar  aaaaar  co  cha  abova,  wa  daaoca  C **  a a  cha  local  acrala  flald 
dacaralaad  bp  aaaualag  alaadc  aggro gaca  raapoaaa  co  aacroacrala  c  aad  dafloa 

“V 

a  fuBCdoo  T  of  pool c loo  within  V  auch  chac 


( 


(a) 


% 


•  T  c 
%  % 


(3.14) 


(2.10) 


c  -  /  ((#>dV  .  /  T  dV  •  I 

%  y  ^  y  %  -V 
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(3.15) 


latroduclag  kiaeaatlcally  idalMlblt,  laflalteslaal  dUplicmat  aad 
correspond  lag  etrala  fields  4uS*c  ,  6(5  due  to  la  carnal  slip  aad  self- 
•trelalag. 


AC  •  ♦  d(S(c)  ,  (3.19) 

aad  a  ralacad  stress  field  c**C\ 

tS(€)  -  t  -  C*1  C(#)  .  (3.17) 

■v  *v  %  \ 

we  find,  fro*  (2.10).  (2.11).  (3.14)  aad  (3.15) 

/  «CS(°dV  •  0  ,  (3.19) 

V  \  -v 


/  tS(c)dV  •  o  -  (/  C“l  T  dV)  .  c  .  (3.19) 

y  \  \  V  '  '  ' 


The  la versa  elastic  coapllaaca  aatrlx  of  the  aggregate  obviously  la 


cro 


C“l  T  dV  . 
%  \ 


Us  also  have,  froa  (3.1)  aad  (3.3), 


(3.20) 


/  «tS(0)dV  •  Ac  -  (/  C  ▼  dV)  .60 

y  \  \  y  '  %  % 


(3.21) 
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•o  chat  Che  aggregate  compliance  matrix  can  ba  altarnately  expressed 


(L  ■  /  C  t  dV  .  (3 

OUcro  x  \ 

Prom  (3.19-3.22)*  (apparent)  macroscopic  plastic  strain  and  "plastic  stress" 
Increments  ere  appropriately  defined  aa 


6cP  -  /  «CS(o)dV  -  6c  -  /  C  dV 

%  y  \  \  y  \  \ 

6oP  -  -  /  «CS(c)dV  -  /  C_1«((i)dV  -  6o  . 

■V  y  X  y  ^  X  X 


(3.23) 

(3.24) 


(For  clarity*  these  various  tsrms  ars  intarprsted  for  the  uniaxial  case  in 
Figure  1.)  Lastly*  we  rslats  ths  Internal  stress  fislds  CS(‘>  «,d  <S(E>  from 

X  X 

(3.3),  (3.14),  (3.17)  and  (3.22): 


;S  (a) 


<S«>  +  C_l  T(c  -  C, 

•C  K  %\  ^Mecro 


o) 


(3.25) 


4.  GENERAL  CHARACTERISTICS  OF  SUBSEQUENT  YIELD  SURFACES 

Hill's  proof  [25]  of  gsnsrallssd  normality  for  aggregates  satisfying  (2.9) 
is  baaed  upon  postulated  conditions  which  are  equivalent  to  the  requirement 

«Tk  6y.  >0  (4.1) 

cr  k 

(l.e.,  crystal  grains  strain  harden  in  activs  systems,  with  icr  denoting  the 
critical  shear  stresa  in  ths  kth  crystallographic  slip  system).  The  final 
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form  of  Hill's  equation  [25]  can  be  written 


jfc  A 

6o  •  6e  -  5e  •  6a  <  0  (4.2) 

f\,  %  % 


In  which  6a,  6e  are  macro-stress  and  related  strain  Increments  which  produce 

^  'V, 

A*  Ac 

slip  In  one  or  more  slip  systems  while  6a  ,  6e  correspond  to  purely  elastic 

'V. 

At  Ac 

response  of  the  aggregate  (with  6e  *  ^acro  ) *  Thus,  we  have  the  following 
orthogonality  conditions  for  the  quantities  defined  in  (3.23)  and  (3.24): 


6a 

'b 


6ep  <  0 


6t 

a. 


6aP  < 

Tj  — 


(4.3) 


(4.4) 


According  to  (4.3),  the  domain  of  all  possible  incremental  plastic  strain  vectors 
6eP  from  a  stress  point  a  in  macrostress  space  is  orthogonal  to  the  domain  of 

%  'b 

* 

all  stress  increments  6a  producing  purely  elastic  response  (i.e.,  directed  into 

'b 

the  elastic  region  in  stress  space).  According  to  (4.4),  the  domain  of  all 
possible  incremental  "plastic  stress"  vectors  6op  from  a  strain  point  e  in 

'b  % 

macrostrain  space  is  orthogonal  to  the  domain  of  all  strain  increments  6c 


producing  purely  elastic  response  (i.e.,  directed  into  the  elastic  region  in 
strain  space).  The  first  of  these  interpretations  is,  of  course,  compatible 
with  the  customary  manner  of  describing  yield  surfaces  and  is  physically  more 

appealing  since  cp  is  the  strain  remaining  upon  mechanical  unloading  to  zero 

% 

stress  (if  the  present  yield  surface  encloses  the  origin  of  macrostress  space). 
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A  particular  yield  hyperplane  In  stress  space  Is  defined  as  (3.1,  3.3  and 


f.  *  N.  I  0  +  -  Tk  -  0 

k  ^  ^  k  cr 


(4.5) 


with  the  corresponding  hyperplane  in  strain  space  given  by  (3.8,  3.14  and  3.17) 


g.  -  N  C-1  T  c  +  t^c)  -  xk  -  0  . 

<\,k  %  /\,  \  k  cr 


(4.6) 


The  shear  stresses  and  are  related  through  (3.8)  and  (3.25): 


Tk<<0  -  +  &  f 1  l  -  «M.cro  «>  •  <‘-7> 


*  * 

Admissible  stress  and  strain  increments  do  ,  £e  corresponding  to  unloudlng  from 
all  potentially  active  systems  must  satisfy 


4  l  1  0 


(4.8) 


N.  C“X  T  5e  <  0 

'uk  \  — 


(4.9) 


which  are  equivalent  statements  since,  from  (3.1),  (3.14)  and  (3.22) 


S  l  •  l  &tac»  • 


(4.10) 
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In  seres*  specs  chs  local  elastic  region  lias  within  the  pyrasld  of  innar 
bounding  hyperplanes  f ^  whose  unit  normal  vectors  era 

b  •  <h  ’ >T/»  \  l  ‘I  •  <*•“» 

The  dlstanca  to  the  ktb  hyparplana  is  (4.5  and  4.11) 

<‘u> 

which  for  Isotropic  sggragstas  slapliflas  to  (sae  Appendix) 

\  •  /T  <tjr  -  t*(c))  .  (4.13) 

Consider  now  the  change  in  position  of  an  active  plane  with  incraaalng  pleatit 
deformation.  Since  k  corresponds  to  the  active  sense  of  slip  in  a  particular 
crystallographic  system,  we  redefine  it  as  (k*)  and  write 

6 Vo  ’  (Vr  ’  6,k(C>>/H?k  X  11  *  (4,U) 


From  (3.13) 


/  ([  «t?(c)  6>.  )dV  «  0  .  (4.15) 

V  k  * 

As  this  Inequality  must  hold  locally  for  the  majority  of  active  systems,  we 
can  assume  it  holds  for  the  particular  k  of  interest.  Thus 
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6t 


S  (o) 
k 


0 


(4.16) 


sod  (4.1  asd  4.16) 

%<♦)  *  4  U^<°>l>/ll?k  l  II  ”  0  *  (A*l7) 

Slallarljr,  th*  disease*  (la  cha  opposite  d  trace  loo  is  serass  apses)  to 
parallel  hyperplaae  corraspoodlsg  to  tb*  negative  Sanaa  of  slip  (k-)  la  this 
crystallographic  spates  is 

"(k->  •  <'«  '  ♦  !k  ’  H  <*•*•> 


and  its  chaags  is 


o 

60 


(k-) 


(4.19) 


Thus,  if  the  cheaps  la  crystal  harden lag  6t 


(k-) 


cr 


la  reverse  slip  la  lass  than 


the  cheaps  la  the  resolved  shear  stress  due  to  latarnal  slip  sad  so lf-s tress lag, 
o 

6d^_j  *  0  sad  the  tvo  hyperplaaes  nova  la  the  sees  diractioa.  This  is 

certaialy  the  case  when  the  crystal  strsla-softeae  la  reverse  slip  (as  la 

suggested  by  the  experimental  vork  of  Pstersoe  (27)  m  copper),  and  there  sill 

be  a  correspond  lag  strong  launch  lager  affact  la  asc  roe  tress  specs.  Thar*  sill 

be  at  least  s  weak  Is  use  h  lager  affact  If  oaly  <.  ,  hones  6D^_)  < 

0 

60 


latroduclag  the  generalised  Schald-Taylor  lav  (24,  28) 


6i  ,  •  1(Y)  6t 

\  X  X 


(4.20) 
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(In  which  H(v)  In  Cb«  crystal  bsrdsnlng  matrix),  this  Inequality  Is  Mt  for  both 
Taylor  hardening  (1,  29)  and  ths  translational  hardening  adopted  by  Budlansky 
and  Uu  (10)  and  Tung  and  Lin  (19) ,  aa  wall  aa  for  any  positive  combination  (30) 

H  •  h(y)l  ♦  c(y)N  NT  .  (4.21) 

(1  Is  an  N  by  N  Mtrlx  all  of  whoso  alsmenta  ara  unity  and  h  and  c  ara  scalar 

\ 

hardening  functions  detaralned  fro*  slngls  crystal  tests.)  In  addition*  H  of 

\ 

(4.21)  Is  at  least  positive  aemi-definlta  (H*0),  and  If  c  |  0,  H>0over  active 

V”  \ 

sysreas  (24,  30).  The  property  M*0  trill  be  of  laportance  In  astabllshing 
uniqueness  of  solution  for  the  discrete  modal  In  Section  8. 

3.  SELECTION  OF  MODEL  FOE  POLY CRYSTALLINE  AGGREGATE  ANALYSIS 

If  we  wars  to  consider  only  (Idealised)  elastically  Isotropic  crystals  in 
defining  an  aggregate  boundary  value  problen,  elastic  field  solutions  for 
point  body  forces  could  be  Introduced  (as  In  (17-18*  20) ),  thus  permitting  both 
non-uniform  microstress  and  displacement  fields  over  unit  cube  faces  A&  while 
still  satisfying  the  virtual  work  aquation  (2.9).  To  consider  aggregate  of 
anisotropic  crystals,  however,  (thereby  enabling  Investigation  of  the  effects 
of  texturing  on  macroscopic  yield  surfaces,  for  exaaple).  It  Is  almost  mandatory 
that  a  modal  with  either  uniform  tractions  or  (at  most)  linearly  varying 
displacements  over  the  fscss  A^  be  adopted.  The  latter  Is  chosen  herein  as 
the  preferred  approximation  on  the  basis  of  the  following  argument. 

Conalder  a  thln-wallsd  tube  subjected  to*  say,  axial  load  and  internal 
pressure.  The  wall  thickness  of  specimens  studied  experimentally  in  combined 
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•treat  tests  Is  often  In  the  range  1-2  non,  with  from  10-30  grains  through  the 
thickness  (see  [31]  and  [32]).  Thus,  as  an  Idealization  of  the  physical 

3 

situation,  we  assume  a  thickness  of  1  mm  and  define  a  unit  cube  V  »  1  mm 
containing  on  the  order  1000  crystal  grains  in  the  corresponding  "flat  sheet" 
representation  (l.e.,  a  macroscopic  plane  stress  problem).  Then  the  longitud¬ 
inal  feces  (Figure  2)  become  planes  of  symmetry  in  our  model  of  Identically 
deforming  cubes.  We  further  assume  the  distribution  of  crystal  orientations 
to  be  symmetric  with  respect  to  transverse  planes.  Hence,  (2.3)  and  (2.4) 
are  setlsfied,  with  the  displacements  either  constant  or  (for  other  combined 
loadings)  linearly  varying  over  the  appropriate  faces.  If  uniform  tractions 
had  been  Imposed  over  the  cube,  the  transverse  and  longitudinal  faces  would  no 
longer  be  planes  of  symmetry  and  adjacent  cubes  could  not  deform  Identically. 
Thus,  we  select  as  a  model  for  analysis  a  unit  cube  (of  generally  anisotropic 
crystals)  on  each  of  whose  faces  either  inf lnlLesimal  displacements  are 
prescribed,  to  give  the  appropriate  macroscopic  strain  increment  (through  2.6), 
or  tractions  ere  zero  (free  face). 

6.  DISCRETIZATION  OF  BOUNDARY  VALUE  PROBLEMS  OF  PRESCRIBED  AGGREGATE 

MACROSTRAIN 

To  discretize  the  above-defined  aggregate  model,  we  Introduce  a 
kinematically  admissible,  approximating  Infinitesimal  displacement  field  which 
Is  continuous  throughout  the  aggregate  and  piecewise  linear  within  each  crystal 
grain.  Correspondingly,  a  crystal  sub-volume  with  constant  microstrain  field 
is  represented  by  a  tetrahedral  element,  herein  called  a  crystallite,  with 
nodal  points  I,  J,  K,  L  (Figure  3).  The  Infinitesimal  displacement  6u(x) 

X  X 

within  the  crystallite  q  is  readily  expressed  in  terms  of  the  nodal 


-M 


displacements  6u  as  (see  [33],  for  example) 


6u(x)  -  l  <J>M(x)  6uM 
^  M(q)  M  ^  * 


(6.1) 


wherein 


>h(«> 


■  “M  +  Vj 


(6.2) 


with  the  general 
in  the  Appendix. 


expressions  for  a^,  in  terms  of  nodal  coordinates  x^  given 
From  (2.2),  the  local  strain  field  is 


M(q)  £q 


M 


(6.3) 


in  which 


8"  ■  f  * 


M(C> 


(6.4) 


Then,  from  (2.10), 


/  6C  dV  -  l  5 l  .  V  -  5c  (6.5) 

y  'v  q  ^(q)  q  'v 

exactly  if  the  infinitesimal  nodal  displacements  on  the  faces  of  the  unit 
cube  are  prescribed  consistent  with  (2.3)  and  (2.6).  Denoting  incremental 
tractions  on  the  crystallite  faces  by  6t  ,  the  virtual  work  equation  for  the 
element  is  (from  2.1) 
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/  6c  .  6£  dV  -  /  6t  •  6u  dS  -  0  (6.6) 

x,  ^  *  e  ^  ^ 


so  that,  from  (6.3)  and  (6.4),  the  element  equilibrium  equation  can  be  expressed 


l 

M(q) 


*1 


dV  -  l  / 

M(q)  Sq 


6t  dS  . 


(6.7) 


Since  an  approximating  displacement  field  has  been  adopted,  tractions  cannot 
be  matched  exactly  between  adjacent  crystals.  Rather,  we  introduce  from  (6.7) 
t’ie  equivalent  nodal  force  Increments 


/  $  T  6t  dS 

S  J  * 

q 


(6.8) 


and  require  that 


I 

q(J) 


0 


(6.9) 


over  the  q  elements  having  the  common  node  J.  Equilibrium  between  crystallites 
is  then  satisfied  in  an  average  sense,  and  the  discontinuous  microstress  field 
6(  through  the  aggregate  is  only  approximately  statically  admissible.  Hence, 
defining  macroscopic  stress  according  to  (2.11),  the  aggregate  virtual  work 
equation  (2.9)  will  not  be  exactly  satisfied,  and  a  measure  of  the 
discretization  error  is  (from  6.5) 


e  •  /  6c  •  6C  dV  -  (/  6;  dV)  .  6e 

y  X  y  f\j 


(6.10) 


which  can  be  written 


It 


*  “  l  6ltai  •  V«  ■  <1  6C/4.v  V  )  .  6c  .  (6.11) 

q  ^iq)  *v(q)  q  jj  %(q)  q  x 

Upon  substitution  of  (6.7)  and  (6.8)  Into  (6.9),  ths  nodal  squllibriua  squat loo 
baconas 

l  <eJ)T  /  it  dV  -  0  .  (6.12) 

q(J)  v"  Vq  ' 

Sines  our  boundary  valua  problsa  Is  ons  of  prsscrlbsd  aacros train,  wa  saparata 
tha  strass  flald  as  In  (3.17)  and  writs  (using  6.3  and  6.4  and  da la ting  tha 
suparscrlpt  (e)  for  simplicity) 


/ 

V 

q 


dV 


M(q) 


flM  «;<•>" 
xq  \ 


V 

q 


«CS  dV 


(6.13) 


wharain  is  tha  crystal  complianca  matrix  rafarred  to  tha  unit  cubs  axas 

and  is  tha  (kinaaatlcally  adnissibla)  lnflnltaslB.il  displacaaant  (laid 

dataminad  by  assuming  alastlc  aggragata  rasponsa  to  6^.  Substituting  (6.13) 
into  (6.12),  tha  nodal  aquillbriua  aquation  is  saparabla  into  two  aquations 
by  daflnltlon  of  tha  flald  6u^*^: 


l 

q(J) 


(fl->T?(J)  I  £  6“(,)H  Va  «  « 

*(q)  N(q)  ^  q 


%q 


(6.14) 


l  /  «?S  <*v  *  ?  •  (6-15) 

q(J)  9  Vq 


Assambling  tha  first  of  thasa  aquations  Into  a  ganaral  matrix  aquillbriua 
aquation  for  tha  ovarall  vactor  6U^*  •  (...,  du* .  ..)*  of  alastlc  nodal 

dlsplacamants  and  noting  that 


a 
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at  c‘l 

*(q)  *c 


(6.16) 


(Mten  £  U  the  cryiul  co^llioci  aatrlx  referred  to  the  crystal  uei  and 

*< 

U  e  transformation  Matrix  given  4 a  the  Appendix),  we  have 

?T?  I  Vc  4?<#>  *  ?  (6.17) 


la  which  ell  crystallites  have  been  choaeo  of  equal  volune  V  «  1  for 
coavealeace  la  subsequent  analysis.  (Thla  la  easily  realised  geometrically 
•lace  e  cubic  voluae  caa  be  separated  Into  alx  equal  volune  tetrahedrons ,  and 
the  unit  cube  can  be  divided  la to  aa  many  sub-cubes  aa  deelred.)  The  Matrix  £ 
la  composed  of  6  by  3  elenenta  1  ,  defined  aa 


0  If  J  la  not  a  node  of  (q)  . 


(6.18) 


The  Matrix  jg  la  a  dlagooal  Matrix  of  positive-definite  aubaatricea 


.-1 


*  •  fi£lJ 


(6.19) 


The  overall  Matrix  l7>  |  la  posit ive-semidef in ite  aloce  a  plecewlae  linear 
dlaplaceaeat  field  admit a  laf lnl tee Inal  rigid  body  notion  with  ao  change  in 
atrala  energy,  la  terne  of  the  uakaown  Interior  diaplaceneata  '  end  the 
preacrlbed  exterior  diaplaceneata  6Jj°  ,  (6.17)  can  be  written 


ll  I  I,  vc  ‘«<#>  *  I  “<#>  -  -  l[  I  «.  *« 
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(6.20) 


wtMrsln  I  it  oktilMd  fro*  B  by  delating  all  vactor  column*  corresponding  to 
exterior  oodti  J°.  Dm  matrix  1  It  similarly  data  reload  by  da la t log  all  vactor 
column*  corraapoodlog  to  lntarlor  aodaa  J.  Dm  aggragata  alaatic  "stlffnass" 
matrix  f  la  aymmatrlc  and  posltlv*-d* finite,  at  arbitrary  rigid  body  motion 
In  tha  solution  having  baao  a  1 lain* tad  by  proscribing  6yc  on  tha  surfaca  of  tha 
unit  cub*.  Dm*,  tha  lntarnal  strain  flald  is  da  tarn  load  as 

•  „{„<  C  •<£’  -  £  il i «. v.  1 

In  uhlch  |C1  danotas  a  row  vactor  of  3  by  3  submstrlcas  of  ^  and  tha  saparata 
contributions  of  lntarlor  and  axtarlor  sodas  of  (q)  ars  as  lndlcatad.  Dm  ovar- 
all  vactor  of  Inf lnltaslmal  stralas  6^**  •  (...,  ,  ...)T  caa  ba  axprassad 

‘I**’  •  iT  n  -  <il  1 1,)'1  ll  fl  I.  *li#  <*•»> 


wTMraln  ^  Is  an  ldaotlty  matrix  and  ^  Is  a  diagonal  matrix  of  alamsnts 
7.  CQfEBAL  SOLUTION  FOB  DfCXDtOfTAL  CRYSTAL  SHEARS 

From  (3.6),  (3.7),  (3.16)  and  (3.17) 

•  h5<v*  «({)  /  <*i‘  -  ?!,» «i» dV  « 

S  R 


or,  using  (6.3), 


* 

) 

* 

1 


l s?s  dV '  f<i>  „L  £  v«  •  «ii>  Kw  l i;  dV 

q  s  q 


,-i 


»“(S)M 


-1  T 

H  *  M* 
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(7.2) 


.5 

in  which  6u  is  the  previously  defined  (Section  3)  infinitesimal  displacement 

V 

— S  T  -S 

field  due  to  internal  slip  and  self-straining  (i.e.,  6£  -  V  6u  ) .  Thus, 

-  X  ^ 

substituting  into  the  second  nodal  equilibrium  equation  (6.15), 


(  #aJ.T  _-l  r  aM  .-(S)M  „  r  .-1  „T  , 

,h>  V  <«»  hI,  ^  4v  V-  *  ,1,  ‘V  f«>  *<’>  i  ^ 


dV  . 


(7.3) 


-fS)M° 

Since,  frosi  (3.18),  6y  0  for  exterior  nodes  on  the  cube  faces,  we  have 

(proceeding  as  before) 


»I  ?  ?,  Vc  *»*  -  ?  ?T  s‘  Vc 


(7.4) 


in  which  6AS  •  (...,  6u^M,  ...  )T,  6y  •  (...,  6y.  ...)T,  and 

x  x  \  x(q) 


?  •  i «(,)«(,» J  *  r 


(7.5) 


where  is  the  transformation  matrix  from  the  local  crystal  axes  to  the 
crystallographic  slip  systems  Hence 


6«s .  <r 1  s  *T «;  »c 


(7.6) 


and 


6C?.  -  I  8"  C  V 

x(q)  M(q)  ™  vi  v  '  c 


.M  -1  _T  „  „T 


(7.7) 
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Upon  subsclcution  of  (7.7)  into  (7.2),  the  internal  stress  field  due  to 
crystal  slip  and  self-straining  becomes 


(7.8) 


-S  -S  T 

Thus,  the  overall  vector  of  incremental  stresses  61  *  (...,  6; ,  .,  ...)  can 

\  'v(Q) 

be  expressed  (from  6.18,  7. A,  7.5  and  7.8) 


61S  ■  -*T  i?  -  ?  ?i  .?  ?i)_1  $1  «T  <7-»> 


with 


(7.10) 


as  previously  defined.  Then,  from  (3.8),  (7.5),  and  (7.9-7.10),  the  overall 

-S  -S 

vector  of  incremental  shear  stresses  6t  -  (...,  6t ,  , ,  ...)  due  to  Internal 

^  ^(97 

slip  is  given  In  terms  of  the  vector  of  infinitesimal  crystal  shears  by 


■  -  *  «  ?T  ■  -  Is 


(7.11) 


wherein 


S  ■  ?  9  -  «i  («i  ?  h*’1  £  *1  •  (7-l2) 


S  T 

The  symmetric  matrix  £  »  ^  Q  ^  will  be  called  the  self-straining  matrix.  In 

Section  8  it  is  proved  that  this  matrix  is  positive-definite  over  critical 
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(potentially  active)  slip  systems,  hence  It  has  a  unique  decomposition  (or 
Inverse) . 

Substituting  (7.9)  and  (6.22)  Into  (6  13)  (or  3.17)  and  denoting  • 

-  X 

(...,  » i  ...)  ,  ws  have 

^(q) 

si  ■  *T  *  u  -  «i<?I « v'X  *v°  -  «T  6i>  (711> 

or  (from  7.12) 


*?,  -  *T  «  l?o  SV°  -  f  •  <7-W) 

The  appropriate  static  admissibility  of  this  stress  field  for  arbitrary  surface 
displacements  and  Internal  shears  Is,  of  course,  confirmed  by  substitution  of 
(7.14)  Into  ths  nodal  equilibrium  aquations  (6.12),  expressed  In  general  matrix 
form,  vhsnce 


4«°  -  «T  sl  i ;  ?  <7-i5> 

T 

since  ^  3  £  from  (7.12).  (The  Influence  matrix  g  Is  analogous  to  the 

lntegrodlfferentlal  operators  on  6£p(jj)  In  [20]  and  [34].)  Ths  overall  vector 
-  -  x 

6t  •(...,  6t  ...... )  of  incremental  rssolvsd  shsar  strsssss  in  the  various 

%  -v(q) 

crystallographic  slip  systems  within  the  cubs  Is  (3.8,  7.5,  7.10  and  7.14) 


•  «  9  l?0  **#  -  «T  sj>  • 


(7.16) 


? 


1 


2} 

Therefor*,  the  equation  for  the  incremental  etrees  in  the  kth  critical  alip 
eye tee  of  cryetelllte  (9)  becoeee 

4w  •  £  C  w«° 

-  i1  ‘i!  V  -  1$  C  (<e)I  4i(,)-  <7-,7> 

For  an  active  ayatea, 

4TM,)  *  4’«,>  *  lk^>4»(,)  <7“> 

with  the  hardeolat  matrix  of  the  cryatailite  depending  only  upon  the  locel 
plaatlc  deformation  (H^  denoting  the  kth  row  vector).  Thu*,  from  (7.17)  and 
(7.18) 


*  <£  £V>T>4i<,)  -  £  C  h,>  <  "i1  %r  ?  vc  ?T  4i 

•  «£  w<°  4<f  -  <  is1  il  i  \ «.  <’  <7 

in  which  the  only  unknown*  are  the  Incremental  aheare  in  the  varloua 

active  alip  aye ten*.  Writing  the  equality  (7.18)  for  all  the  active  oyatema, 
equation*  (7.19)  can  be  expreaaed  in  general  matrix  fore  a* 


'«*  ♦  <*  «*  "J1 4u  *  i!*  «*  <«.»*  4«° 


(7.20) 


a 


tfhmlB  ^  li  a  dlaiooal  aatrli  vtwia  tlmoti  ait  tba  individual  cryatalllta 
hardaslat  nstrlcss.  Iba  subscript  A  dsaotss  ths  rsspsctlv*  vsctor  or  aatrla 
rsducad  to  lacluds  only  tboss  crystallltss  (q)  contslnlnf  00s  or  aor*  sctlvs 
critical  systsas.  Tbs  synastrlc  aatrli 

^  ? A  *A  (7*2l) 

Is  posltlvs-daflalta  (or  sll  crystal  bar dsn Inf  astrlcss  of  typ«  (4.21) ,  ss  pro vs a 
la  ths  follow Inf  ssctlon.  sad  (7.20)  (or  7. If)  ylslds  s  unlqus  solutloo  for  4v  . 


I.  UflQUEXESS  Of  IfTEKXAL  FIELDS  If  THE  D1SCKET1ZED  ACCtECATE  MODEL 

Froa  (6.14),  (6. If)  snd  (7.10),  ths  overall  vsctor  of  laf laltsslasl  slsstlc 

•a  .  T 

atralaa  dE  ■  (»..,  •••)  Is 

*  £  I*  $  *  <••» 

■slatroduclaf  ths  scalar  svsrafs  u#  of  (3.f ) ,  wo  writs 

=.  •  i  •*<,>  •  S(,)H(,>  ve  •  >i  ■  *r v,  *  °  •  <,j> 

Hsacs,  froa  (7.14)  sad  (8.1), 

Z9  -  (|o  4¥°  ’  ?T|J>  Vc>  0  .  (8.3) 

AaauBs  bow  two  distinct  sots  of  latsraal  f  Is  Ids  *X(q) 
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both  °*  **lch  MClafy  all  appropriate  equations  of 
tho  discretized  aodel  for  <o  lof lnlceslaal  aacros train  6c  (l.o.,  uniquely 
prescribed  6^°).  Dsootlng  tbolr  differences  by<6£^p»«  6t^j  -  6(^j, 
etc.,  then  fro*  the  ebove 


H(q)*  4^(q) * 


<6t >7  N  $  $T<6^>  >  0  .  (8.6) 

This  equation  is  written  only  over  potentially  active  slip  sy stone  corresponding 

to  the  current  scstes  of  Internal  stress  and  strain,  since  the  active  systeas 

for  either  set  of  lncrenental  shears  will  belong  to  this  critical  group. 

S  T 

Therefore,  f  -  ^  Q  N  Is  positive  definite  over  critical  systeas;  consequently. 
It  Is  also  positive  definite  over  active  systeas  (J£$  >0).  We  also  find, 
froa  (3.7),  (7.5),  and  (7.16), 

I <««,,)> •  ♦  «j>T  mii1  <*i> ■  ° •  <»•*> 

Introducing  an  Inequality  due  to  Hill  (24]  written  over  critical  systeas  of  s 
crystal  grain. 


<tlw>  •  <4S(,>>  i  «(,>  <4i(,)>  •  «•*> 


and  substituting  Into  (8.S),  we  have 

0>  r^(q)J<«Y>  ♦  <«^^^T<6y>.  (8.7) 
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Fro*  (8.4)  the  second  tern  Is  positive,  and  for  jg  >_  0  (as  In  4.21)  the  first 
term  Is  non-negative.  Hence,  the  Inequality  Is  violated  for  non-zero 

'V. 

from  which  we  conclude  thst  the  lncrenental  sheers  ere  unique.  It  follows 
from  (7.14)  that  the  Incremental  stress  field  Is  unique.  Consequently,  the 
Infinitesimal  total  strain  field  Is  elso  unique,  and  the  proof  Is  complete. 

The  matrix  of  (7.21)  Is  obviously  positive-definite  from  (8.4)  and  the 
above. 

A  final  comment  Is  nscessary  psrcaining  to  determination  of  the  unique 
aggregate  response  from  e  perticular  deformed  state,  since  an  admissible 
solution  of  (7.20)  Is  constrained  by  the  physical  requirement  6y^  >_  0  for  all 
k  (which  Is  Implicit  In  8.6).  If,  for  example,  (7.20)  Is  solved  for  6y  based 

'X/ 

upon  the  expectation  that  all  critical  slip  systems  will  be  active  In  a 

macros train  Increment  6c  ,  certain  of  the  6y.  ,  »  may  be  calculated  as 

'v  xiq; 

negative.  In  this  event  (7.20)  must  again  be  solved,  after  eliminating  the 
appropriate  slip  systems,  until  all  Incremental  plastic  shears  are  positive. 

9.  REMARKS  ON  QUANTITATIVE  MODELS  AND  THE  CALCULATION  OF  AGGREGATE 
YIELD  SURFACES 

At  a  particular  stage  of  aggregate  straining,  a  subsequent  yield  surface 
Is  obtained  by  determining  the  positions  of  ylsld  hyperplanes  In  stress 

T 

space,  as  follows.  Consider  the  case  of  applied  biaxial  strain  c  -  (e. .,€__) 

^  XX  *A 

-  -  -  x 

corresponding  ‘:o  the  biaxial  macrostress  stats  o  ■  (ain,o__)  discussed  In 

\  ii  a 


Section  5.  From  the  elastic  solution  of  (6.21)  (or  6.22)  we  define,  according 
to  (3.14), 


(9.1) 


iV 
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In  which  j  Is  a  6  by  2  matrix  whose  columns  are  determined  from  separate 
evaluations  of  (6.21)  for  the  states  (c^2.*c22  *  ^  end  (e(  .  ■  0,  e^7)  respec¬ 
tively.  Hence,  from  (3.20) 


I  c~fL\  T,  N 

q  ^(q)  ^(q) 


V 

c 


(9.2) 


a  2  by  2  matrix  (using  only  the  first  two  rows  of  each  j£^j).  Adopting  the 
volume  average  definition  of  macroscopic  stress  (2.11),  we  have 


o 

'V. 


ioi) 


(9.3) 


where  again  only  the  two  leading  elements  are  chosen.  Introducing  the 
notation  cR  for  (since  from  (3.1)  and  (3.3)  iB  the  residual  micro- 

X  X  X 

stress  field  remaining  in  the  aggregate  upon  unloading  to  zero  macrostress), 
then,  from  (4.7),  (6.16)  and  (7.5), 


-R  -S  c  -I  -  -  . 

Tk(q)  Tk(q)  +  ^k  *?c  ^(q)  ^  ^Macro  £  ^ 


(9.4) 


-S  -s 

wherein  t.  .  .  is  the  sum  of  increments  6t,  .  .  determined  from  (7.11)  over  the 
k(q)  k(q) 

strain  history  and  C„  and  o  are  as  calculated  above.  The  distance  to  the 

'V'Nacro  x 

kth  hyperplane  in  macrostress  space  then  is  determined  from  (4.12)  as 


Dk(q)  "  (T«q>  "  Tk(q))/H  ^k(q)  *(q)  H 


(9.5) 


in  which,  from  (4.10),  (6.16)  and  (7.5) 


>ili 


* 


? 


-1 


^k(q)  ^(q)  ^k  £c  £(q)  £(q)  ^Mai 
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(9.6) 


The  direction,  as  given  by  (4.11),  remains  fixed.  The  elastic  domain  Is  the 
Inner  bound  of  all  such  hyperplanes  in  a^,  space. 

A  quantitative  model  which  should  prove  suitable  for  such  evaluations 
(corresponding  to  thin-walled  tubes)  is  an  aggregate  cube  of  1296  crystallites 
arranged  within  216  sub-cubic  volumes.  From  the  symmetry  conditions  of 
Section  5,  the  discrete  boundary  value  problem  can  be  reduced  to  a  consideration 
of  324  crystallites  with  112  distinct  interior  and  exterior  nodes  and  223  unknown 
nodal  displacements.  The  symmetric,  positive-definite  matrix  is  then  of  order 
223  with  a  half-band  width  of  approximately  43.  Thus,  its  accurate  inversion  is 
a  relatively  modest  task  for  a  triangular  decomposition  routine.  Concerning 
the  solution  of  (7.20)  for  incremental  shears,  the  order  of  should  be  small 

(less  than  100)  during  the  early  stages  of  aggregate  straining  as  only  a 
relatively  few  slip  systems  will  be  critical.  With  increasing  plastic  deforma¬ 
tion,  numerical  studies  could  feasibly  be  continued  until  several  systems 
become  active  within  each  of  the  majority  of  the  crystallites, at  which  stage 
the  order  of  £A  would  still  be  less  than  1000.  A  digital  computer  program 
for  studying  the  quantitative  effects  of  aggregate  texturing,  crystal  structure 
and  anisotropy,  and  crystal  hardening  laws  on  theoretical  macroscopic  response 
is  in  preparation,  with  numerical  results  to  be  presented  in  subsequent 
papers. 
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APPENDIX 

Internal  stress  and  Infinitesimal  strain  tensors  are  expressed  In  vector 

form  as 


?  m  ^  2  (^2*  ^ ^22*  ^  ^  ^23*  ^33^  (A.l) 

65  -  (65n,  /T  6512,  /T  6513,  6522,  /T  6523»  6533)T  (A.2) 

from  which  the  equilibrium  and  kinematic  aquations  can  be  written  V  £  ■  £ 
and  PX65  ■  6ji  =  (6u^,  6u2»  6u^)X,  with 
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0  0 


i  »j  0  0  0 

0  >2  ^J3  0 

<iJi  0  *3 


(A.  3) 


(wherein  ^  denotes  partial  differentiation  with  reepect  to  the  corresponding 
spatial  coordinate). 

The  kth  row  vector  of  the  crystal  transformation  matrix  £  Is  given 

k  k 

terms  of  unit  vectors  ^  ,  A  In  the  normal  and  gilds  directions,  respectively, 
of  the  kth  crystallographic  slip  system: 


"k 


(♦Jxj  +  *2X1  ^  (#ix3  +  *3*1  )  * 

♦***.  +  ♦jlj  )  . 


(A.4) 


Consequently 


(A.  5) 


For 

V  Vl 

summing 


the  tetrahedral  crystallite  (Figure  3), the  configuration  parameters 
are  determined  as  (Introducing  the  permutation  tensor  end 
on  repeated  Indices) 


a 
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1  J  K  L 

6Vq  eijk  Xi  XJ  *k 


1  K  I  L 

6V"  Cijk  Xi  Xj  *k 


1  I  J  L 

6V  eijk  Xi  Xj  *k 


1  I  K  L 

6V  €ijk  Xi  XJ  *k 


(A.  6) 


and 


®Mj  "  6V"  <cof*ctor>M 


j 


1 

1 

1 

1 


whnra  6V  la  tha  datarnlnant  of  tha  matrix. 


I 

*1 

J 

*1 

K 


(A.  7) 


Tha  traniformatlon  matrix  relating  the  atreaa  or  lnflnltealmal  atraln 

components  referred  to  the  unit  cube  axes  with  those  referred  to  the  crystal 


*(q)  " 


axes  (Figure 

3)  la 

2 

*11 

n  *11*12 

^  *11*13 

2 

*12 

r/2  a12a13 

2 

*13 

* *11*21 

*11*22+*12*21 

*11*23+*13*21 

^  *12*22 

*12*23+*13*22 

^  *13*23 

^  *11*31 

*11*32 +*12*31 

*11*33+*13*31 

^  *12*32 

*12*33+*13*32 

*  *13*33 

- 

2 

*21 

*  *21*22 

^  *21*23 

2 

*22 

^2  *22*22 

2 

*23 

^  *21*31 

*21*32+*22*31 

*21*33+*23*31 

SI  *22*32 

*22*33+*23*32 

*  *23*33 

2 

*31 

S2  •31a32 

^  *31*33 

2 

*32 

^  *32*33 

2 

*33 

(A.  8) 


- awtaz+zx*  . 


•  • 


In  which 


COB  (X^ 


Thu* 


$<q)  "  £(q)$(q)  ’  6£<q)  "  £(q>  6£<q> 


35 


(A.  9) 


Incremental  Macroscopic  Variables 
for  tha  Uniaxial  Casa 


Longitudinal  plants 
of  symmatry  ^ 


Fig.  2 


Idtallxtd  Unit  Cuba  in  Thin-Walltd 
Polycrystalllnt  Sptcimtn 
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Fig.  3  Tstrohsdral  Crystallite  Showing  Oritntation 
of  Crystal  Axss  and  Slip  Systom 
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